Phonon-assisted relaxation kinetics of statistically-degenerate excitons 

in high-quality quantum wells 
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Acoustic-phonon-assisted thermalization kinetics of excitons in quantum wells (QWs) is developed 
for small concentrations of particles, p2u J; 10^ cm~^, when particle-particle interaction can be 
neglected while Bose-Einstein statistics already strongly influences the relaxation processes at low 
temperatures. In this case thermalization of QW excitons occurs through nonequilibrium states and 
is given by the following scenario. During the first transient stage, which lasts a few characteristic 
scattering times, the correlations with an initial distribution of QW excitons disappear. The next, 
adiabatic stage of thermalization usually takes many characteristic scattering times, depends only 
upon two control parameters, the lattice temperature Tt, and the degeneracy temperature Tq oc p2D, 
and is characterized by a quasi-equilibrated distribution of high-energy QW excitons with effective 
temperature T{t). We show that the thermalization law of high-energy particles is given by ST{t) = 
T{t) - Ti, oc e~^°^/t, where Ao is a marginal value of the continuous eigenvalue spectrum of the 
linearized kinetics. By analyzing the linearized phonon-assisted kinetics of statistically-degenerate 
QW excitons, we study the dependence Ao = Ao(Ti,,To). Our numerical estimates refer to high 
quality GaAs and ZnSe QWs. Finally, we propose a special design of GaAs-based microcavities, 
which considerably weakens the bottleneck effect in relaxation of excitons (polaritons) and allows 
us to optimize the acoustic-phonon-assisted thermalization processes. 

PACS numbers: 78.66-w, 72.10.Di, 63.20.Kr 



I. INTRODUCTION 

The formation, resonant or phonon-assisted, of QW 
excitons and their following relaxation towards a fi- 
nal (quasi-) equilibrium thermodynamic state at low 
lattice temperature T;, are |±he._sU|bject of numerous 
experimental and theoreticaBa'DLl'lju studies. Recently, 
relaxation thermodynamics has been formulated and de- 
veloped in order to analyze how Bose-Einstein statis- 
tics of high-density QW excitons mfluences the phonon- 
assisted thermalization processes.l3 The above thermo- 
dynamics assumes a strong, dominant exciton-exciton 
scattering and, therefore, relaxation through quasi- 
equilibrium thermodynamic states. The relaxation ther- 
modynamics has successfully been applied to model the 
P2D-dependent thermalization and photoluminescence 
kinetics observed in early experimentsBllj with high- 
density excitons (5 x 10^ cm~^ < p2D < 10"'^^ cm~^) 
in GaAs QWs.eI 

The long-lived indirect excitons in high-quality 
GaAs/AlGaAs coupled QWs provide a unique oppor- 
tunity for studying quantum degeneracy in a system of 
two-dimensional bosons. In this case, the long radiative 
lifetimes of indirect excitons allows the system to cool 
down to temperatures where ^hr .'T^'ll^^^fi' M^M^^r" S^-s be- 
comes statistically degenerate.yyy'BBnO The qual- 
ity of present day GaAs/AlGaAs coupled QWs has been 
considerably improved in conpnarisnn ..witb. those used 
in the pioneering experimentsBBlHjEjB'EHl one decade 
ago. Furthermore, the very recent magneto-optical 
experimentf£j clearly indicate that the in-plane momen- 
tum ?ik|| of indirect excitons is a well-defined quantum 
number in high-quality GaAs/AlGaAs coupled QWs. 

Thermalization of hot photoexcited excitons down 



to the temperature of the cold lattice occurs mainly 
via scattering by thermal bulk longitudinal acoustic 
(LA) phonons and is much more efficient for quasi-two- 
dimensional (quasi-2D) systems as compared to bulk 
semiconductors. This follows from the relaxation of mo- 
mentum conservation in the z-direction (the QW growth 
direction) for quasi-2D systems: the ground-state mode 
k|| = 0, i.e., the energy state E = 0, couples to the 
continuum energy states E > Eq, rather than to the 
single energy state E ^ Eq — 2Mxv1 {vg is the longi- 
tudinal sound velocity and is the in-plane transla- 
tional mass of excitons) as occurs in bulk materials. As 
a result, the LA-phonon assisted kinetics of QW excitons 
becomes dominant at p2D ^ 1 — 3 x 10^ cm~^: in this 
case exciton-exciton scattering can be neglected while 
Bose-Einstein (BE) statistics already strongW influences 
the relaxation process at low temperatures^ Crossover 
from classical to quantum statistics occurs near the de- 
generacy temperature fcsTo — iTrh^ P2d/ {gM^), where 
g is the spin degeneracy factor. For p2D = 3 x 10^ 
cm^^ the degeneracy temperature of indirect excitons 
in GaAs/AlGaAs coupled QWs is Tq = 0.79 K. This 
estimate refers io g = 1, which can be achieved in 
the 6-type GaAs / AlGaAs jCoupled QWs by applying a 
static magn|etic field H|jz.B Note that the very recent 
experimental^ deal with GaAs/ AlGaAs coupled QWs at 
extremely low cryostat temperature Ti, — 0.05 K. 

In this paper we study analytically and model nu- 
merically the acoustic-phonon-assisted relaxation kinet- 
ics of statistically-degenerate QWj_excitons at low den- 
sities. The recent experimentaHfl'E3 allow us to visual- 
ize, by means of LO-phonon-assisted emission, the LA- 
phonon-assisted kinetics of quasi-2D excitons in ZnSe 
QWs and, in particular, to prove that for p2D ^ 10^ cm~^ 
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the above kinetics indeed occurs through non-eqpilihrium 
distributions of QW excitons. In experirnentsEjOo'tj 
the concentration of BE-degenerate indirect excitons in 
GaAs/AlGaAs coupled QWs usually varies from p2D ^ 
10^" cni~^ at the very end of an optical excitation pulse 
to (02D ^ lO'* cm^^ at large delay times t ^ 50 ns after 
the optical excitation. Thus, thermalization of the in- 
direct excitons at large delay times cannot be described 



within the relaxation thermodynamicda and does need a 
separate theoretical analysis. 

The classical Boltzmann kinetic equation has been 
generalized in order to iaclude quantum statistics by 
Uehling and Uhlenbeck.Es The relevant quantum ki- 
netic equation for spatially-homogeneous dilute system 
of statistically-degenerate quasi-2D excitons coupled to 
bulk LA-phonons is 



d 



q 



+ [A^k||<(l + ^k||+qJ-(l + iV^)(l+nf)iVk||+q„]5(i?k|, -i?k,,+q,, 



(1.1) 



r 



where A^ky and are the occupation numbers of ex- 
citon in-plane mode ky and phonon bulk mode q — 
{q||,9z}, respectively, and qy is the in-plane projection 
of q. The terms in the first and second squa re brack- 
ets on the right-hand-side (r.h.s.) of Eq. (1.1) describe 
the Stokes and anti-Stokes LA-phonon-assisted scatter- 
ing processes, respectively. The acoustical phonons are 
assumed to he in a thermal equilibrium at the bath tem- 
perature TfjJj The matrix element is given by M{q, q^) — 
[{Dlhq)/{2pv,V)]^/'^F^{q^Lj2), where p is the crystal 
mass density, Dr^ is the deformation potential of exci- 
ton - LA-phonon interaction, Lz is the thickness of a 
QW, and V is the volume. The form- factor Fz{x) = 
[sin(x) /x] [e*-^/(l — x^/7r^Uj-efers to an infinite rectangu- 
lar confinement potential.EZl The latter function describes 
the relaxation of the momentum conservation law in the 
z-direction and characterizes a spectral band of bulk LA- 
phonons, which effectively interact with QW excitons. 
Note that Eq. (^]^) is valid only for the kinetic stage of 
thermalization, Lep-bfifore a low-temperature collective 
state of excitonsEEaHy builds up. 

The main aim of our work is to study the fundamental 
features of the acoustic phonon-assisted thermalization 
kinetics of QW excitons from initial strongly nonequilib- 
rium A^ky (t = 0) towards the final equilibrium distribu- 
tion with well-developed Bose-Einstein statistics, when 
n ^ 1- Our numerical simulations of the LA- 

k||=0 ~ ^ 

phonon-assisted kinetics clearly demonstrate that after 
the first transient, which lasts a few characteristic scat- 
tering times, a slow adiabatic stage of thermalization 
builds up (see Fig. |l]). This stage is characterized by 
a quasi-equilibrium distribution of high-energy QW exci- 
tons with effective temperature T{t) = Ti, + ST{t) and is 
independent of the initial distribution at t = 0. The adi- 
abatic stage lasts many characteristic scattering times 
and arises due to the need to populate the low-energy 
in-plane modes with A^k||~o ^ 1 in the presence of ef- 
fective suppression of the stimulated scattering processes 
(e.g., an intense incoming. Stokes flux of excitons into the 



ground-state mode ky = is nearly compensated by the 
outgoing, anti-Stokes scattering out of the state ky=0). 
In order to describe the adiabatic stage of relaxation, we 
find a generic solution of the quantum kinetic Eq. ( |1 . l| ) . 
While the generic solution we calculate is different from 
that derived in Ref. for the phonon-assisted kinetics of 
bulk excitons at Tf, < Tc (Tc is the critical temperature 
for Bose-Einstein condensation of excitons in bulk semi- 
conductors), similar to this case it depends only on two 
control parameters of the system, Tq and T^. Further- 
more, the only gross information from a particular shape 
of the initial distribution at f = is absorbed by the start 
time tc of the adiabatic stage of thermalization. This is 
shown in the Fig. |l|(a), where we vary the parameters of 
the initial Gaussian distribution. In turn, the first tran- 
sient depends upon initial distribution and lasts only a 
few scattering times Tsc, as illustrated in Fig. |l|(b). 

Because at any nonzero bath temperature > the 
occupation number of the ground-state mode ^k||=o ~ 
exp(ro/Tb) — 1 is finite, the adiabatic stage ends up with 
linearized kinetics. The latter kinetics refers to the last, 
exponential stage of thermalization at t — > cxd and is char- 
acterized by A^k||=o(t) - ^k||=o T{t) - Tfc cx e^^"*. 
For statistically-degenerate QW excitons the thermaliza- 
tion time Tth = l/Ao(T{,,ro) reaches its smallest values 
at bath temperatures E^/kB ^ Tfc ^ Tq. While the 
above inequality does not usually hold in the experi- 
ments with statistically-degeBei[at p |i ndir | eC| t excitons in 
GaAs/AlGaAs coupled QWsJia-yyyyilll we describe 
a special design of GaAs-based microcavities (MCs) for 
optimization of the LA-phonon-assisted thermalization 
kinetics of low-density MC polaritons. In the proposed 
microcavities with a large positive detuning between the 
cavity and QW exciton modes, the MC polaritons have 
radiative lifetimes on a 100 ps - 1 ns time scale, so that the 
phonon-assisted relaxation towards well-developed Bose- 
Einstein statistics with large occupation numbers can op- 
tically be visualized. Thus, the MC design we discuss is 
an interesting alternative to the semiconductor micro- 
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FIG. 1; (a) Transient relaxation towards the adiabatic stage 
of evolution calculated for various initial Gaussian distribu- 
tions N^{t = 0) (X e-i5 '525(E:-e)^ ^ijgj.g £ ^ ?i2kjj/(2Af^£;o), 
and parameter £ = 3 (dash-dotted lines), 4 (solid lines), and 
5 (dashed lines). The gross dependence of evolution upon 
initial conditions is absorbed by time ti ~ tc- (b) The first 
stage of LA-phonon-assisted relaxation {t < Tsc) for a partic- 
ular distribution of indirect excitons with e = 4. In both plots 
Tf, — 0.25 and To = 1.14. The dimensionless values, time in 
Tac and energy/temperature in Eq, can easily be rescaled to 
the dimensional units by using Eq — 33 neV and Tsc = 41 
ns for GaAs/AlGaAs coupled QWs, and Eq = 162 /leV and 
Tsc = 2.7 ns relevant to ZnSe QWs. 



cavities with zero detuning between the cavity mode and 
QW excitons, where huge nonclassical occupation num- 
bers of the lowggg/-^^ggolariton states have recently 

In numerical evaluations we use = 0.21 mo, 

Vs = 3.7 X 10^ cm/s, and = 15.5 eV, relevant 
to GaAs/AlGaAs coupled QWs, and = 0.86too, 
Vs = 4.1 X 10^ cm/s, and — 6.9 eV, relevant to single 
ZnSe QWs, respectively (toq is the free electron mass). 
Note that the disorder-induced scattering and localiza- 
tion processes, not included in our model, are relatively 
strong in up-tcudate ZnSe-based QWs and require a seft 
arate analysisHH However, in the very recent worktJ 
the first fabrication of high-quality MgS/ZnSe/MgS QWs 
with less than 1 ML fluctuations of the well width and, 
therefore, with extremely low inhomogeneous broadening 
has been reported. 

In Sec. the Boltzmann equation is adapted in or- 
der to formulate the acoustic-phonon-assisted kinetics of 
statistically-degenerate QW excitons. We also discuss 
some approximations for the form- factor Fz{x)i which 
describes the relaxation of momentum conservation in 
scattering of QW excitons by bulk LA-phonons. 



In Sec. Ill, we find the generic solution of the acoustic- 



phonon-assisted kinetics from a strongly nonequilibrium 
initial distribution of QW excitons Ng{t = 0) towards a 
final Bose-Einstein distribution with large occupation 
numbers of the low-energy states, N^^^ > 1, where e is 
the dimensionless energy defined by e = h\y{2M^Eo). 
The generic solution is independent of a particular shape 
of the initial distribution and describes the adiabatic 
stage of thcrmalization, which starts at t = tc in a few 
characteristic scattering times after t = 0. We show that 
at i 3> tc thermalization of high-energy QW excitons is 
given by ST{t) cx e~^°*/t, where Xq is the lowest positive 
eigenvalue of the relevant linearized kinetics. 

In Sec. the linearized phonon-assisted kinetics of 
statistically-degenerate QW excitons is formulated and 
analyzed. We show that the eigenvalues {A} of the lin- 
ear collision integral form a continuous spectrum, oo > 
A > Ao > 0, separated from the nondegenerate eigen- 
value A = 0, and that the corresponding eigenfunctions 
{'(/'e(A)} have three well-defined isolated critical points. 
The dependence of the marginal eigenvalue Ao on the 
bath (Tb) and degeneracy (To) temperatures is studied. 

In Sec. straightforward numerical simulations of the 
phonon-assisted relaxation kinetics at Tf, ^ Tq are com- 
pared with the generic solution of the quantum Boltz- 
mann equation. We also show that at the beginning of 
the adiabatic stage of thermalization, at tc < t ^ , 
the population dynamics of the ground-state mode is 
given by N^^oit) oc (l+xt)'^ , where the parameters x and 
ly are calculated analytically. Furthermore, we propose 
a particular design of GaAs-based microcavities (MCs), 
which compromises the efficiency of LA-phonon-assisted 
scattering (the density of states is cx M^) with the degen- 
eracy temperature To cx M~^, i.e., allows us to avoid the 
bottleneck effect in relaxation and, therefore, to optimize 
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the thermalization kinetics of low-density QW excitons. 

In the Appendix, some relationships relevant to 
the thermalization dynamics of quasi-equilibrated high- 
energy QW excitons are given. 



II. BOLTZMANN KINETIC EQUATION FOR 
DEGENERATE QW EXCITONS 

For hot QW excitons, which are in-plane isotropically 
distributed a,t t — 0, the thermalization kinetics due to 
bulk LA-phonons can be treated in one-dimensional en- 
ergy space [see Eq. (2) of Ref. In the following we 
express energy E and temperature T (as well as Tt, and 



To) in terms of Eq, i.e., we use the dimensionless values 
oi E ^ e = E/Eq, and T {kBT)/Eo. In order to 
derive and analyze the generic solution for relaxation at 
Tb < To I it is convenient to rewrite the above equation 
in terms of the variable 



fe{t) 



(2.1) 



where Ns{t) and iV° = - 1] are the current 

and (final) equilibrium distribution functions of QW ex- 
citons, respectively, and the chemical potential fi is given 
hy II — Ti)ln{l — e~'^°^'^''). In this case the kinetic equa- 
tion reduces to 



at Tsc 



/•Ss(e) /-co 

/ Fs(e,ei)£s(e,ei,*)£idei + / -FAs(e, ei)'CAs(£, ^i, O^i^^ei 



(2.2) 



where the Stokes (S) and anti-Stokes (AS) collision inte- 
grands are 

Cs{e,e,,t) = [/,(t)~A_,,(t)](l + nPh + 7vO_^J 

+ T,(iVt,J7e(i)/e-s,(t), (2.3a) 

£As(£,£l,t) = [A(t)-/e+e,(i)]«:'-iV,VJ 

-r,(A^O+,J'/e(t)/s+e,(0- (2.3b) 

The distribution of thermal bulk phonons is given by the 
Planck formula, i.e., n?'' = l/(e^/"^'' — 1). The scattering 
time is defined by t^c — {n^U p)/{DlM^Vs)- The func- 



tions 9g /AS (^) I which determine the integration limits on 
the right-hand-side of Eq. (2.2), are 




0, s< 1/4 

2^^- 1, 1/4 < e < 1 

e, 1 < e, 

1-2^, e< 1/4 

0, e>l/4. 



(2.4a) 
(2.4b) 



The functions -Fs/as(£:£i) i^^ Eq. (2^) are given by 



„S/AS 



-Ps/As(£:£i) 



,S/AS 



\F;:{aeia)\" 



-1/2 



da, 



(2.5) 



where u] 



S/AS 



[1 - (V^ - V^TiT) Ve?] and d, 



,S/AS 



1 /2 

[l — {^/e + ^/F^Ti)"^ / ef\ . The dimensionless param- 
eter a in the argument of the form-factor function is 
given by a = {LzMxVs)/h. 

In order to estimate the contribution of the form-factor 
Fzix) to t he s pectral functions Fg/ji^g{e, ei) one can an- 
alyze Eq. ( ^.5|) for e — > 0. In this case only the spe ctra l 
function Fas(£,£i) is relevant to the kinetic Eq. (2.2), 
and Eq. dJ) yields 



F., 



a\/ei{ei - 1) 



(2.6) 



For ei > 1, the spectral width Ae (FWHM) of Fas (0, ^i) 
is given by Ae ~ 2.26/a and determined solely by 
Fz{a£i)- For thickness = 8 nm of GaAs QWs (see 
Ref. |l^) one gets a ~ 0.054 and, therefore, Ae ~ 42. 
The latter estimate can be rewritten in the dimensional 
energy units as AE = EqAe = A.52fiVs/Lz ~ 1.44 meV. 
Thus the spectral band of LA-phonons, which scatter 
a low-energy QW exciton, can be evaluated as AE ^ 
VsApz, where Apz ~ h/Lz is the uncertainty of the mo- 
mentum in the z-direction due to the QW spatial confine- 
ment. The above estimates show that for relatively cold 
QW excitons with energies e < Ae (the effective temper- 
ature T''^ = AE/ks ^ 16.5 K for Lz = 8 nm) the form- 
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factor can be approximated by Fz{aeia) = Fjfi) = 1- 
In this case the integral on the r.h.s. of Eq. (2.5) can be 
calculated explicitly: 



Fs/As{£,£i) - ^sTas 



(2.7) 



where 



m = 



z {f [arcsin (VC) , t] -^(f)}' 0<e< 1 



(2.8) 



Here, F{(f>, m) = [l - msin2(6l)] dO and K{m) 
F{-K/2,m) are the elliptic and the complete elliptic inte- 
grals of the first kind, respectively. 



III. GENERIC SOLUTION 



{5T=^,N=Nl] 




= 0, Nit = 0)} 



In this Section the thermalization kinetics at Tf, < 
To is described in terms of a generic solution, which 
weakly correlates with the initial distribution (i = 0) . 
In particular, we derive a thermalization law for high- 
energy excitons and examine nonequilibrium distribution 
of low-energy QW excitons. We also specify a reference 
point for the generic solution, i.e., a set of parameters, 
which unambiguously determines the calculated evolu- 
tion. Schematic picture of the relaxation kinetics in phase 
space is shown in Fig. ||. The generic solution is relevant 
to the times t > tc, where tc is the start time of the 
adiabatic stage of evolution. 

As was emphasized in Sec. I, the generic solution as- 
sumes a quasi-equilibrium distribution of high- energy ex- 
citons {e > 1/4): iV,>i/4(i) = l/[e(^-m)/Tit) _ -^i^^ ^j^jg 

quasi-equilibrium distribution is characterized by the ef- 
fective time-dependent temperature T{t) — Tb + ST{t) 
and chemical potential fl(t) whose time variations are 
supposed to be small, so that //(i) can be taken the same 
as the chemical potential of final equilibrium distribution 
p,{t) ~ II. Indeed, for T{t) < Tq the chemical potential 
is given by fl{t) ~ -T{t)e^^°/^'^^'^ < T{t) and, therefore, 
= \il{t) — ii\ <^ 5T{t). Since at the final stage of re- 
laxation kinetics the effective temperature of high-energy 
excitons approaches the bath temperature {ST{t) — + 0), 
starting from some moment in time one meets the con- 
dition 5T{t) < Tb (provided that > 0). Therefore we 
can linearize the quasi-equilibrium distribution function 
with respect to ST{t). In this case one obtains 



L>i/4{t) = -[ST{t)/T^]e. 



(3.1) 



On the other hand, as we show below, a low-energy 
kernel (0 < e < 1/4) of the distribution function charac- 
terizes the non-equilibrium QW excitons. With increas- 
ing time t > tc, the low-energy kernel shrinks in energy 
space, and the ratio |A'^j<i/4 — -^°<i/4l/^"<i/4 decreases 
approaching the limit when the linearized kinetics be- 
comes valid. 



FIG. 2: Schematic picture of the relaxation kinetics at 
Tb ^ To. The reference (start) point {to, STo, Nq^} unambigu- 
ously determines the calculated (solid line) evolution from the 
initial noequilibrium distribution Ne{t — 0). The dashed line 
corresponds to the generic solution, tc is the start point of 
adiabatic stage of evolution. For the times t — tc > Tsc both 
evolution spirals almost coincide showing a unique path to- 
wards the final equilibrium distribution {5T — 0, Ni: — N^}. 



A. Nonequilibrium distribution of low-energy QW 
excitons 

In order to analyze the evolution of low-energy QW 
excitons within the scenario described in the previous 
subsection, we substitute Eq. ( BT ) into Eq. ( |2.2| ) and get 
the reduced kinetic equation for e < 1/4: 



§-^fe{t) = - Me) + Us)STit)] fS) + v{e)5T{t), 



(3.2) 



where 



4 

eo(e) = — / Fas(£,£i) (n^ - iV!'+,J eidei, 

Tsc Jl-2y^ 



(3.3) 



a(£) = ^ / FAs(£,£i)(iV^+ej'£i(£ + £i)rfei, 

(3.4) 



T]{e) 



TAs(£,£l)«-iVe"+eJ 



xei(e -f ei)d£i. 



(3.5) 
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and Fas{s + £i) is given by Eq. (2/?). Equation (3.2), 
which describes the population dynamics of the low- 
energy states during the adiabatic stage of thermalization 
{t > tc) , is a, linear inhomogeneous differential equation 
for /e<i/4. Its complete solution can be written as a sum 
of the homogeneous and inhomogeneous contributions: 



x6T{T)dT, (3.6) 



where 



fe<T/Ait) - e-(*-*«)««(-)-''(*-*«)«^(-)/,<i/4(i = to), (3.7) 



ST{T)dT, 



(3.8) 



and to is an arbitrary reference (start) time for the cal- 
culated evolution, i.e., to > (see Fig. |^). 

Thus the adiabatic stage of the relaxation kinetics 
into the lower-energy states is completely determined by 
N,{t) = +Tb{N°Yf,{t), provided that one knows the 
reference (start) distribution fe<i/A{t = ^o) f^nd ther- 
malization law 5T — 5T{t) for high-energy QW exci- 
tons. Note that in a sharp contrast with the acoustic- 
phonon-assisted relaxation kinetics at T;, < Tc in three- 
dimensional systemSjCHI the thermalization dynamics of 
low-energy QW excitons depends upo n the homogeneous 
contribution /]l2"y4 given by Eq. (3.7). 

As will be shown in the next subsection, the thermal- 
ization law for quasi-equilibrium high-energy QW exci- 
tons (e > 1/4) is given by 



5T{t) = 



STo 



Al — Ao 



t~tn ' 



(3.9) 



where Aq = ^o(O) characterizes the inverse thermalization 
time att~^ oo, Ai ~ ^o(l/4) is another characteristic pa- 
rameter relevant to the beginning of the adiabatic stage 
(Al :s> Aq), and STq = ST{t = to) determines reference 
(start) effective temperature T{t — tp ) = Tt + STg for the 
calculated evolution. Equation (3.9) is valid for t > to. 



Using the thermaliz atio n law (3.9) we find the integral 
on the r.h.s. of Eq. (jsj): 



p(t,ii) 



STo 



^ {Ei[-Xo{t-to)]- Ei[-Xoiti~to)] 
Al — Ao 

-Ei[-\i{t - to)] + Ei[^Xi{ti - to)]} , (3.10) 



J^(e */t)dt is the exponential integral 



where Ei{z) 
function. 

The reference distribution fe<i/4(t = to), which deter- 
mines the homogeneous solution ( |3.7[ ), is independent of 
the initial distribution of hot QW excitons, Ne{t — 0). 
There is, however, an integral relationship between STo 
and ,fg<i/4(t = to), which we will discuss in subsection 
[lie along with possible approximations for f^<i/i{t — 
to)- 



/e<i/4(i) = /e<V4(i) + vi^) / e-f*-^)?"*^)-''**^^'^^^^) tion dl) 



B. Thermalization of high-energy QW excitons 

The temperature change ST ~ 6T{t) characterizes the 
time evolution of high-energy quasi-equilibrated particles 
through Eqs. (|2.lD and ([3.1]). In order to derive the tem- 
peratur e la w (B.gi) we substitute fe<i/4 and /e>i/4 given 
by Eq. (3.6) and (3.1), respectively, into the kinetic equa- 



dt 
where 

lit) 



ST{t) 



[ao + ai{t)]ST{t)+pST^{t) 



4Tl 



1/4 



Fs{ec,ec-e) (l -t 
xf^{t){ec - e)de, 



-lit). 

(3.11) 
(3.12) 



and the parameters ao, (3, and functional ai (t) - 
Q^i [/e< 1/4(0] defined in Appendix. Equations ( 3.11 ) 



( |3.12D refer to some energy 1 from the high-energy 
domain e > 1/4. 

In the adiabatic stage of relaxation, when t > tc, one 
has ST{t)/Tb <C 1 so t hat a i and PST^ can be neglected 
on the r.h.s. of Eq. ( 3.11 ), because |Q;i(i)| <C ao and 



\PST^{t)\ < \lit)\. At the end of the adiabatic stage 
it — tc ^ Ao^^) the distribution function Ni.{t) is already 
very close to 7V° even for small energies e < 1/4 and, 
therefore, the phonon-assisted relaxation kinetics be- 
comes exponential, i.e., ST{t) cx e~^"'. For statistically- 
degenerate QW excitons (iV°^Q 1) one estimates from 
Eqs. (3.3) and (Al) that positive Ao = C (0) is much less 
than positive ao. In this case Eq. ( ^.11 ) can be solved 
iteratively. The first iteration, which can formally be ob- 
taine d by neglecting the time derivative on the r.h.s. of 
Eq. (pl|), yields 



ST{t) = 7(0/«o. 



(3.13) 



For the same time domain i — tc ^ Aq ^ Eq. ( |3.6| ) yields 
the following approximation: 

h<i/Ait) - e-(*-*o)«o(-)/at - to) + ^,ST{t). (3.14) 

40 (e) 



By substituting Eq. (|3.14|) into the r.h.s. of Eq. ( |3.12| ) 
one derives from Eq. pll): 



1 /■i/'* 

STit) = ^J Fsisc, Sc -e){l + ?.P^,Af°) /e(to) 
xe-(*-*o)«o(-')(£c-e)de, (3.15) 

where the constant a is defined in Appendix. We can 
further simplify Eq. ( 3.15 ) taking into account that ^o(e) 
is a monotonously increasing function of energy and that 
at large t only small vicinity of e = contributes to the 
integral. Finally, we end up with the asymptotic law: 



7 



1 0.0205 
1 0.02 

\ , M 0.0195 

\ ^ 
'\ 0.019 

'\ 0.0185 




1 

\ = 0.25 






100 200 300 400 500 





0.005 
0.004 
0.003 
0.002 
0.001 


100 200 300 400 500 

FIG. 3: Time dependence of S T an d A^^^i (inset) calculated 
by using the thermalization law (3.9) (dashed lines) and by di- 
rect numerical modelling of the initial kinetic Eq. (2.2) (solid 
linesV The control parameters Tb and To are the same as in 
Fig. |. 



where the constant 70 is giv en in Appendix. Equation 
( |3.16|) is identical to Eq. provided that the reference 
(start) temperature Tj, + ST{t = to) of the high-energy 
quasi-equilibrated QW excitons is determined by 



C. Reference point for the generic solution 



In order to determine the relaxation kinetics of low- 



energy particles, i.e., 



e<l/4 



e<l/4 



tc), one 

needs to specify the reference distribution f^<i/i{t — to) 
and STo = ST{t = to) = T{t = to) ~ T^. " According 
to Eq. ( ^.7| ), with increasing time t — to only the small 
vicinity of e = gives contribution to the homogeneous 
part of the generic solution. Therefore we approximate 



£<l/4 



(t = to) 



Ng^oit = to). The above ap- 



STo = 6T{t = to) = (7o/5)(Ai - Ao). 



(3.17) 



proximation assumes that the reference time to > tc is 
close enough to the start time tc of the adiabatic stage 
of relaxation so that 7Ve(io) < iV° at e < 1/4. Thus, 
using Eq. (2.1) we determine the reference distribution 
at t — to by 

7VA" - N° 

/e<l/4(i-to) = ^^. (3.18) 

The latter expression is completely defined by the only 
one unknown parameter Nq"^, the population of the 
ground-state mode at the reference time to- 

Equation (|3.18D allows us to find STo = ST{t = to) 
through the integral relationship ( |3.15 ) taken at t = to- 
Furthermore, within the approximations used in the 
derivation of the asymptotic law ( |3.16 ) a simplifie d (al - 
gebr aic) form of this relationship, given by Eqs. (3.17) 
and (|A5|), is valid as well. After the value of STo is de- 
termined, using the thermalization law (3.9) one can eas- 
i ly fi nd the time dependence of the distribution function 
( |3.l[ ) of high-energy excitons. Therefore, the three pa- 
rameters to, Nq^, and STo, completely specify the refer- 



I I I I ^-^-^ ^ ^\.J J U 7 — — — ^ \J ■} — ^j- J "XT J ^ 

While the above derivation of Eqs. (^JJ)-(|JJ) as- gnce point for the generic solution, as illustrated in Fig. |, 



sumes that t — tc > t ~ to Xn i we have checked nu- 



merically that the temperature law ( p.9[ ) holds through 
the whole adiabatic stage, i.e., for the time interval < 
t — to < 00. Furthermore, the numerical evaluations also 
clearly indicate that the start temperature ST(t = to) is 
practically independent o f the local energy £c ^ 1 used 
in the derivation of Eqs. ( ^.ll| )-( |3?T7| ). In Fig. | we plot 
ST = ST{t) and A^^^i = iVe=i(i) calculated by using the 
thermalization law (|]^) (dashed lines) and b y dir ect nu- 
merical modelling of the initial kinetic Eq. ( |2.2|) (solid 
lines), respectively. 

For the time domain X^^ < t — to Xq^ the thermal- 



ization law ( |3.9D can be approximated by ST{t) = (Ai — 
Xo)~'^STo/{t-to). The latter dependence ST{t) cx l/{t- 
to) is consistent with that found for the adiabatic stage of 
the phonon-assisted relaxation kinetics of 3D bosons (ex- 
cittps) at Tb <Tc, when Bose- Einstein condensate builds 
up.E2l For 2D systems at nonzero the occupation num- 
ber of the ground-state mode N^J^o is always final. There- 
fore the exponential kinetics ST{t) oc e~^°^*~'^°\ which 
results from Eq. (3J) for t — to^ Xq^, develops at the fi- 
nal stage of relaxation, when |A^e=o(i) — -^£=ol/^£=o ^ 1- 
Note that because for statistically-degenerate QW exci- 
tons, when Tb < To and N^^q 1, one has Ap ^ ^ Ts, 



IV. LINEARIZED KINETICS FOR 
STATISTICALLY DEGENERATE QW EXCITONS 

If for any energy e the distribution function N^^ of 
quantum-degenerate quasi-2D excitons is close enough 
to final so that \N^ - N^\/N^ < 1 and fs{t) becomes 
small, the phonon-assisted kinetics can be linearized. In 
particular, the adiabatic stage of thermalization at times 
t — to > Aq ^ refers to the linearized kinetics. The lin- 
earized kinetics can be described in terms of the real 
eigenvalues {A} (A > 0) and the corresponding eigen- 
functions {^£(A)} so that /^(t) = CAi/'e(A) exp(— At). 
The initial kinetic Eq. ( |2.2| ) reduces to the linear Fred- 
holm equation of the second kind with respect to '0e(A): 



A7Ae(A) = — 



where 



es(e) 



eAs(e) 



Fs(e,ei)£s(£,ei)£irfei 



FAs(e, ei)'CAs (e, £i)£ic!ei 



(4.1) 



(see Section IV), the two modes of behavior, ST on 1/t 
and ST cx e~", are well-separated in time. 



£s(£,ei) = mX) ~ V's~e,(A)] (1 + nf^ + N^_,J , 

(4.2a) 
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FIG. 4: A typical shape of the eigenfunction The 
main part of the figure shows the first order pole, which arises 
at the energy region e < 1/4 (A — 0.0043/rsc in this particular 
example). The top inset illustrates eigenfunction behavior at 
critical points in the high energy band e > 1. The eigenvalue 
spectrum (A in units I/tsc) is shown in the bottom inset. The 
control parameters Tt and To are the same as in Fig. hi. 



£As(e,ei) = [i^eW - V-e+ei (A)] {n- 



ph 



(4.2b) 



and functions and 0as are given by Eqs. (2.4a) and 
( 2.4b| ), respectively. 



Thus we replace the solution of Eq. (2.2) by the eigen- 
function analysis of the Fredho lm i ntegral equation (4.1). 
The numerical solution of Eq. (4T ) clearly shows that for 
a given Tf, > all except one eigenvalues {A} are non- 
degenerate, positive, and belong to the continuous spec- 
trum. This is illustrated in the bottom inset of Fig. ^, 
where the set of eigenvalues shown by the stars covers 
the same interval cxd > A > Aq = Ao(Tb,ro) more dense 
the more discrete points in energy space are used. The 
isolated non-degenerate eigenvalue A = is due to conser- 
vation of the total number of QW excitonsJn our model 
(the only integral of motion of the systemo) . 

Since all {A} are nondegenerate, the corresponding 
eigenf unctions {'(/'^(A)} form a basis in Hilbert energy- 
spac e. F or the energy band < e < 1/4 one derives from 
Eq. (|t]): 



£<l/4 



(A) 



a{e) 



A-eo(£)' 

where cr(e) is a smooth regular function of e given by 



(4.3) 



aie) 



'sc Jo 



Fs(e,ei)V'c+ei(A) 



and ^o(e) is defined by Eq. ( |3.3[ ). Thus the eigenfunc- 
tion tpeW has an isolated singularity (first order pole) 
at ex — C(7^(A). The singularity is integrable in terms of 
principal value integration. A typical shape of the eigen- 
functions at e < 1/4 is shown in Fig. |[ In the energy 
band e > 1/4 the eigenfunction has another singularity 



FIG. 5: The inverse thermalization time Ao (in units l/rs< 
as a function of the control parameters Tt and To. 



at the point (1 — y^) . This singularity is logarith- 
mic, i.e., integrable. It arises when the singularity of 
'0e-ei(A) at e — ei — e\ [see Eq. (4.2)] coincides with 
the upper boundary of integration, 0s{e), in the Stokes 
collision term on the r.h.s. of Eq. (4.1). The step-like 
jump at the critical point (1 -j- ^Je\] 



^does not accom- 
panied by discontinuity of the eigenfunction. The latter 
critical point originates from a singular behavior of the 
form-factor function Fs{e,ei) at £i = 2y/e — 1. The fea- 
tures of the eigenfunction V'e(A) at points e — {Izk^/ex)'^ 
are shown in the top inset of Fig. ^. 

The marginal point Ao = Ao(Tb,ro) of the continu- 
ous spectrum of {A} is indeed given by Aq = fo(0): for 
A — > Ao the singularity point e\ 0, i.e., approaches its 
lowest limit. Because A^^^ is the longest relaxation time 
generated by the continuum Aq < A < oo, the eigenavalue 
Aq determines the phonon-assisted kinetics at < — s- oo 
and, in particular, yields the characteristic thermaliza- 
tion time in the relaxation thermodynamics developed in 
Ref. ^. The dependence of Aq upon the control parame- 
ters of the system, T}, and Tq, is plotted in the Fig. ^. 

The eigenvalue Aq can reach both limits, i.e., Ao S> 
and Ao <C (see Fig. |5|). The first limit, which corre- 
sponds to the considerable acceleration of the thermal- 
ization kinetics in comparison with that in the 3D case, is 
due to the relaxation of momentum conservation in QW 
exciton - bulk acoustic phonon scattering. The slowing 
down of thermalization, Aq «C t~^, occurs (i) at low bath 
temperatures T;, < 1, due to the exponentially decreasing 
number of thermal acoustic phonons with energy e > 1, 
and/or (ii) for well-developed quantum statistics (To ^ 
Tf,), due to an effective suppression of the stimulated ki- 
netic processes. There is no direct phonon-mediated in- 
teraction of low-energy QW excitons £ < 1/4, and at low 
temperatures Tf, < 1 the relaxation kinetics occurs by the 
two-step process: "low-energy QW exciton (e <C 1/4) 
+ phonon (ei ~ 1) ^ QW exciton (e2 = £ + ei ~ 1) 
low-energy QW exciton (£4 = £ + £1 — £3 <C 1/4) 
+ phonon (£3 ~ 1)". The first, anti-Stokes transition 
quenches with decreasing temperature Tf, < 1 and yields 
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a temperature dependent bottleneck effect in thermal- 
ization. In turn, the critical slowing down of the relax- 
ation kinetics at Tq » Th arises due to mutual compen- 
sation of two stimulated fluxes, into and out the low- 
energy QW states e <C 1. For example, for the ground- 
state mode £ = the collision integrand responsible for 



the stimulated kinetics is given by A^e=o(^ei>i 



i)- 



At Tq ^ Tf, the latter combination becomes small in 
spite of a large occupation number iV^^o ^ Ij because 



>1 



\i>ll 



>1 



^£l>ll 



as a result of a 



very small value of the chemical potential, \fj,\ <C 1. 

Now w e can give an alternative prove of the tempera- 
ture law ( |3.9|) which refers to high-energy QW excitons 
with £ > 1/4. Namely, the function f^{t = to) can be 
expanded over the basis {■^'^{X)}. Then the solution of 
Eq. (2.2) is given by 



(4.5) 



where c\ are the expansion coefficients. At large times 
i — io 3> Tth only a small vicinity of A near Aq con- 
tributes to the integral, due to the time exponent in the 
integrand. If now one assumes a regular distribution of 
quasi-2D excitons at t — Iq, the coefficients cx smoothly 
depend upon A. As a result, the approximation c\ ~ c\ „ 
can be used in the integrand on the r.h.s. 



of Eq. (|^ 

Furthermore, the eigenfunctions V'e(A) have nearly the 
same smooth shape at high energies, as illustrated by 
Fig. H Thus we can also put V'e(A) — '0e(Ao) on the 
r.h.s. of Eq. (^). As a result, both cx and V'e(A) can be 
extracted out the integral. Using Eq. (3.1) we immedi- 
ately get ST{t) cx e"^"'*"*") /{t-to), which coincides with 
Eq. ( |3.9| ) at i — to ^ Aj^^. Note that the above derivation 
is based on the particular spectrum (A = -I- contin- 
uum Aq < A < oo) of the linearized collision integral and 
has no analogy in the relaxation kinetics due to particle- 
particle interaction. In the latter case the fivefold degen- 
erate eigenvalue A = is separated from the contiguous 
spectrum by a set of discrete isolated eigenvalues. ell 



DISCUSSION 



In order to test the generic solution (3.3)-( 3T0| ) we 
model the phonon-assist ed r elaxation of excitons within 
the initial kinetic Eq. (^]^) reduced to energy space 
[Eqs. (2.1)-(p^)]. An adaptive inhomogeneous grid with 
100 - 200 points for e is used to cover the close vicinity 
of the ground-state mode £ = (the maximum val ue o f 
the dimensionless energy is e^ax — 20). Equation ( |2.2| ) 
is evaluated by a fourth-order Runge-Kutta integration 
routine with the time step 0.001 — O.OItsc. In order to 



calculate integrals on the r.h.s. of Eq. (2^) we perform 
spline-approximation for fe{t) at every iterative step. 

In numerical simulations we use the dimensionless tem- 
peratures Tf, and Tq and measure time in Tgc- This 
makes our results suitable for various QWs and sets of 



the control parameters, provided that Eq, and Tsc are 
specified. In Fig. ^ time evolution of the distribution 
Ne{t)^ N° + TbiN^Yfeit) as a numerical solution of 
Eq. ( p.2[) is compared with the corresponding generic so- 
"T^) relevant to < Tq. 



lution (|3.3D-( p.lC| ) relevant to < Tq. AU plots demon- 
strate an excellent agreement between analytical and nu- 
merical solutions. Note, that at high Tf,, e.g., Tb — 10 (see 
the top plot in Fig. ^), the thermalization time Tth = Aq ^ 
achieves the limit Tth ^ 'Hsc, and the relative duration 
of the adiabatic stage estimated in terms of tc (dura- 
tion of the first transient) becomes smaller than that at 
Tfc < 1. In this case the influence of the initial distribu- 
tion {t — 0) slightly affects the calculated evolution at 
the beginning of the adiabatic stage, and the analytical 
solution fits numerically simulated data a little worse as 
can be seen for the distribution functions at t = 0.015rsc. 
In contrast, at low the relaxation kinetics at the adi- 
abatic stage is slow. For example, at Tf, < 0.25 it lasts 
more than 1000 scattering times, where typical values of 
Tsc in GaAs/AlGaAs coupled QWs are on a scale of tens 
nanoseconds. 

In Fig. 1^ we also compare time dependences for the 
ground-state mode population N^^q = N^^oit) obtained 
numerically (solid lines) and analytically (dashed lines). 
Aga. in , this figure shows that at e = the generic solution 
(|3.3|) -( PT0|) reproduces the adiabatic stage of the phonon- 
assisted relaxation kinetics very well. Within the time 
interval < t — to < Tth the generic solution yields the 
following simple approximation for the adiabatic kinetics 
into the ground state mode: 



(5.1) 



Parameters x an d v can b e fo und comparing series ex- 
pansions of Eq. (|5.l[ ) and ( ^ ) about the point t ^ tg. 
In such a way we get X = (c? — C2Nq'^)/ {ciNq'^) and 
1/ = cf/(cf — C2N^), where time- independent constants 
ci and C2 are given by 



ci 

C2 



(A^,%-7V^")[Ao 



Wo77(0)(7V,%)', 



1 



x{2A^ + <5Toa(0)[5Ao + Ai + 2(5To6(0)]} 
-r6(5To77(0)(Ar,%)'[3Ao + Ai + 25To^i{0)]). (5.2) 



Time dependences of the ground-state mode popul ation 
Ns=o{t) calculated at = 1 and Tb = 2 using Eq. ( [5T| ) 
are shown in Fig. |^(a) with dash-dotted lines. At t — to > 
Tth the approximation (5.1) violates because the phonon- 
assisted relaxation kinetics star ts to become exponential 
As we have shown in Sec 



< 



109 



IV 



for a given concentration 
the thermalization ki- 



of QW excitons p2D 
netics slows down with decreasing Tb/To, i.e., with devel- 
opment of quantum statistics. By increasing both tem- 
peratures. To and Tf,, and keeping unchanged the ratio 
To/Tb ^ 1 one can simultaneously avoid the above bot- 
tleneck effect in relaxation and achieve high population of 
the ground-state mode, A^e=o ^ 1. However, in this case 
the concentration of excitons p2D oc Tq increases as well. 
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FIG. 6: Evolution of the distribution function A'^^ at the 
adiabatic stage calculated for various sets of the control pa- 
rameters, Tb a nd T q. Solid lines correspond to numerical eval- 
uation of Eq. ( |2.2[ ) ; das hed lines are obtained using the generic 
solution given by Eq. (3.6). 
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FIG. 7: Population dynamics of the ground state mode 
Ne=o ~ Ne=o{t) calculated for different bath temperatures 
at To = 4,10 (a), and To = 1.14 (b). Similarly to Fig. |, 
the solid and dashed curves are calculated by using Eqs. (fiM) 
and (B.6), respectively. Dash-dotted lines correspond to time 
dependences N^^oit) calculated using approximate Eq. ( |5.l| ). 



so that exciton-exciton interaction eventually becomes 
the main mechanism of relaxation in GaAs / AlGaAs cou- 
pled or ZnSe single QWs. 

In contrast, in high-quality GaAs-based microcavities 
with a relatively large positive detuning S = h{uJo — Ut) 
between the cavity mode (luoo) and QW exciton line 
(hiOt) the LA-phonon-assisted kinetics remains dominant 
even at relatively high degeneracy temperatures Tq ^ 
Tf, > 1. In these systems, by changing the detuning S 
within the band (5 = 50 — 100 meV one can design an 
effective in-plane mass Mx so that the bottleneck effect 
in phonon-assisted scattering, due to the low density of 
states cx M^, is already relaxed, whereas the degeneracy 
temperature Tq oc P2d/ (EqMx) is still relatively high. 
Indeed, in such MCs the lower polariton branch gives 
rise to the in-plane translational mass much smaller than 
the mass of optically-undressed QW excitons. In the 
meantime the excitonic component fx of the microcav- 
ity polaritons is already very high, (px ^ 0.999, resulting 
in the long optical decay (in the z-direction) lifetimes 
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FIG. 8: Possible design of GaAs-based microcavities: (a) ex- 
citonic component (fix = i^x(k||) of tlie MC polariton eigen- 
state; and (b) the lower branch polariton dispersion h,{u) — ujt). 
Detuning 5 = 50 meV (dashed lines), 75 meV (solid lines), 
and 100 meV (dash-dotted lines). The energy of ground-state 
QW excitons is given by -Bk||=o = ^(^t = 1.522 eV. 



on a few nanosecond time scale. This is illustrated in 
Fig. ||(a), where the detuning 5 is equal to 50, 75, and 100 
meV. The corresponding polariton (exciton) masses are 
given by — 0.023mo, 0.050mo, and O.OSItoo, respec- 
tively. The relevant lower-branch polariton dispersions 
are plotted in Fig. ||(b). Because the energy i?o is only 
on a 0.01 meV energy scale, the parabolic approximation 
of the lower-branch dispersion curves is valid for low- 
temperature relaxation kinetics. Thus the LA-phonon- 
assisted thermalization of the exciton-like MC pol ar iton s 
can indeed be modelled by the kinetic Eqs. (2.1)-( ^^ . 
Time evolutions of the exciton distribution (t) , which 
are typical for the proposed design of GaAs-based MCs, 
are shown in Fig. ^ (see the plots with Tb — 1, and 10). 
Figure ^(a) illustrates the corresponding population dy- 
namics of the ground-state mode, N^^q = N^=o{t). 



the thermalization kinetics from an initial distribution of 
QW excitons at t = occurs through the nonequilibrium 
states. The following conclusions summarize our results. 

(i) For the case of well-developed Bose-Einstein statis- 
tics, when Tf, < Tq so that N^^q > 1, the relaxation 
kinetics of QW excitons coupled to thermal bulk acous- 
tic phonons is given by the following scheme. Within a 
few characteristic scattering times the correlation of the 
distribution function N^{t) with the initial N^{t = 0) 
disappears, and the subsequent thermalization of QW 
excitons is described in terms of the adiabatic stage of 
relaxation. The adiabatic stage is characterized by the 
start time tc, which absorbs a gross information about 
the initial distribution N^{t — 0), and by the parameter 
Ao, which depends only upon the bath and degeneracy 
temperatures, and Tq. At the beginning of the adia- 
batic stage, i.e., for the time domain < t — tc ^ Aq ^, 



-^°<l/4l 



■^£<i/4' ^^"^ ^^"^ thcrmaliza- 



one has |iVe<i/4- 

tion kinetics is strongly nonexponential, with ST on 1/t 
and A^£=o oc (1 -t- X'^Y ■ At large times, when the de- 
viation of the system from the final equilibrium state 
is already smaU (|iVe(i) - N^\/N^ < 1), the adiabatic 
stage of the phonon-assisted thermalization becomes ex- 
ponential, 5T oc e~^°*, and can be described within the 
linearized kinetic equation. 

(ii) The linearized LA-phonon-assisted kinetics of QW 
excit ons is formulated in terms of the Fredholm integral 
Eq. (4.1). The eigenvalues {A} of the collision integral 
are given by the continuous spectrum Aq < A < oo and 
the isolated eigenvalue A = 0. The marginal eigenvalue 
Ao = Ao(T'{,,To) determines the thermalization time at 
t — > cx) by Tth = ■ In dependence on the two control 
parameters, Tf, and Tq oc p2D, the thermalization time 
achieves two limits: Tth ^ '''sc and Tth 3> t^. The eigen- 
functions ipsi^) of the collision integral (4.1) are smooth 
integrable functions with three isolated critical points in 
energy space. The critical points of tpeW give rise to a 
first order pole, a logarithmic singularity, and a continu- 
ous step-like jump. 

(iii) Because the LA-phonon-assisted kinetics becomes 
dominant only at small concentrations of QW excitons, 
P2D ^ 10^ cm^^, nonclassical statistics of quasi-2D ex- 
citons in ZnSe or GaAs QWs develops at very low bath 
temperatures Tf, < 1 K. The proposed design of GaAs- 
based microcavities with a relatively large positive de- 
tuning between the cavity mode and QW exciton line, 
h{ujo — LOt) > 50 meV, allows us, however, to build up 
-^e~o ^ 1 by means of QW exciton - bulk LA phonon 
scattering in much more favorable conditions, i.e., T;, ^ 1 
K and Tth < Tsc- 



VI. CONCLUSIONS 



In this paper we have studied thermalization kinetics of 
statistically degenerate QW excitons coupled to thermal 
bulk acoustic phonons. For concentrations of QW ex- 
citons P2D ^ 10^ cm~^ the particle-particle interaction 
in GaAs or ZnSe QWs can be neglected in comparison 
with QW exciton - bulk acoustic phonon scattering, and 
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APPENDIX A: TEMPERATURE LAW 

In Appendix we give the exp ressi on s for the parame- 
ters and functions used in Eqs. ( 3.11 ), ( [3.15| ), and ( [3.16 ). 
The parameters ap, and (3 arise from those terms in the 
colhsion integral which contain /e>i/4(t) and /^>i/4(i), 
respectively. Collecting such terms together one gets 



ao 



Be I Fsiec, e) (1 + nf + NI_^) ede 

ec-l/4 



ec-1/4 



POO 

/ FAs{ec,e){nf~Nl+^)e^de 
Jo 



(Al) 



and 



TscTh 



ec-1/4 



Fs(ec,e)(iV°_J'(£e-e)erfe 



FAs(ec,e)«+e)'(£c + e)erfe 



(A2) 



The function ai{t) stems from the terms proportional to 

fe>i/'i{t)fe<i/4{'t) and is given by 

«iW = — Fsie,,e,- sKN^y mis, - s)de. 

Tsc Jo 

(A3) 



The parameter ce from Eq. (3.13) is obtained by col- 
lecting all the terms oc ST{t): 



V{£c - e) 



ede 



Fs{e^,e){l + nf + Nl_^) 
1/4 l;o[£c-£) 

(aoT,c£c) / (4T2) . (A4) 



When deriving Eq. ( p.l6| ) from Eq. (" 3.15) we first put 
e = everywhere in the integrand (3.15) except the expo- 
nent. By changing the integration variable £ to ^ = £ole) 
and putting e = we derive the asymptotic Eq. (3.16), 
valid in the limit t —^ oo. In this equation 70 is a time- 
independent pre-factor given by 



70 



^ Fsiec, ec) (1 



\ £cf £=o{to) 



NLo 



(A5) 
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